Copyright © 1973 American Telephone and Telegraph Company 

The Bell System Technical Journal 

Vol. 52, No. 6, July-August, 1973 

Printed in U.S.A. 



Coupled Mode Theory of Round 
Optical Fibers 

By D. MARCUSE 

(Manuscript received January 12, 1973) 

This paper presents a comprehensive theory of mode coupling in 
optical fibers with imperfections. The paper begins with the derivation of 
a general coupled wave theory based on the modes of the ideal fiber. The 
general theory is applied to a simplified description of guided and radia- 
tion modes of the fiber that is valid for small core-cladding index differ- 
ences. The simplified theory results in expressions for the coupling coeffi- 
cients that are nearly as simple as those of the slab waveguide. As an 
example, the theory is applied to the calculation of radiation losses caused 
by pure core diameter changes and by elliptical deformations of the fiber 
core. 

I. INTRODUCTION 

Dielectric optical waveguides support a finite number of guided 
modes and an infinite number of radiation modes. ' Even if the number 
of guided modes that can be supported by the waveguide is reduced to 
one, the presence of the infinite number of radiation modes forces us 
to be concerned about mode conversion phenomena. Coupling among 
the guided modes of a multimode optical waveguide (multimode wave- 
guide refers to the guided modes) is caused by imperfections in the 
refractive index distribution or the geometry of the optical waveguide. 
Its effects may be beneficial for reducing the delay distortion that re- 
sults from uncoupled multimode operation. 23 Coupling between guided 
modes and the continuum of radiation modes is usually not desired 
unless the waveguide is intended to serve as an antenna. However, a 
certain amount of coupling is unavoidable and results in scattering 
losses. 1 - 4 ' 5 

E. G. Rawson has calculated light scattering from fiber waveguides 
with irregular core surfaces by an approximate technique. 6 However, 
his method is not suitable to calculate coupling between guide modes. 

817 
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A theory of radiation losses in round optical fibers has been presented 
in Ref. 7. This theory was based on a field expansion in terms of the 
exact modes of the fiber. Because of the complicated mode fields, a 
theory based on the exact modes of the guide is very tedious and re- 
sults in equations whose numerical evaluation is difficult and costly. 
However, Snyder 8 and Gloge 9 have shown that the description of the 
modes of dielectric waveguides can be greatly simplified if it is assumed 
that the difference of the refractive indices of core and cladding is only 
very slight. This assumption allows an approximate treatment of the 
mode problem, resulting in very much simpler field expressions. The 
coupled mode theory based on approximate modes yields expressions 
for the coupling coefficients that are almost as simple as those for the 
slab waveguide. The entire coupled mode theory thus becomes simpler 
and its numerical evaluation becomes cheaper. 

This paper starts out with a derivation of the coupled wave equations 
in terms of modes of the ideal guide. This mode description is somewhat 
different from the coupled mode theory in terms of local normal modes 
used by Snyder. 10>n It results in simpler expressions for the coupling 
coefficients. The exact coupled mode theory is then applied to the 
problem of coupling between the simplified waveguide modes. We limit 
the discussion to coupling caused by changes or imperfections in the 
waveguide geometry. Coupling caused by refractive index inhomo- 
geneities, which could be handled in a similar fashion with the use of 
the exact expressions for the coupling coefficients and the approximate 
mode description, is not discussed in this paper. 

Finally, we apply our results to the problem of scattering losses of 
guided modes caused by diameter changes and elliptical deformations 
of the waveguide core. We also derive simplified expressions for the 
coupling coefficients between guided modes far from cutoff and discuss 
coupling between guided modes caused by deformations of the fiber 
core and by curvature of the waveguide axis. 

II. EXACT COUPLED MODE THEORY 

The dielectric optical waveguide with imperfections is defined by 
a certain refractive index distribution n = n(x, y, z) that enters 
Maxwell's equations: 

V X H = *w€ n a E (1) 

V X E = -«a»/ioH. (2) 

E and H are the electric and magnetic field vectors of a general field 
distribution in the waveguide. The fields are assumed to have the time 
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dependence e iut , with radian frequency a>; e and \i are the vacuum 
values of the electric permittivity and magnetic permeability. In addi- 
tion to the refractive index distribution n of the real waveguide, we 
consider the index distribution n„ that defines an ideal guide from which 
the real guide deviates in some way. 

We decompose the fields into their transverse and longitudinal parts. 
The electric field is thus represented by the equation 

E = E, + E z (3) 

and the magnetic field is given by 

H = H ( + H z . (4) 

By using a similar decomposition of the V operator (e* is a unit vector 
in the ^-direction), 

V = Vt + e 'fz' (5) 

Maxwell's equations can be written in the form 

V, X H, + (e 2 X 25A = zW n 2 E< (6) 

V,XE : + (e, X ^) = -zuvzoH,. (7) 



and 



The longitudinal field components are expressed in terms of the trans- 
verse field components 

E z - — ^ : V, X H, (8) 



and 



H, = - J- V f X E t . (9) 

lUHo 



The modes of the ideal waveguide with index distribution n (x, y) are 
defined as solutions of the equations 

V t X M„ - i/3,(e, X SCO = iue n z „S,t (10) 

and 

V t X &,. - ip,(e, X 8,0 = —iuiioXpt. (11) 

The index v is a mode label and /3„ is the propagation constant of the 
j4h mode. The longitudinal components of the mode fields are similarly 
expressed as 

E„ = ^—5 V, X K vt (12) 
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and 

5C„ = - 4- V, X 8*. (13) 

The transverse field of the waveguide with the index distribution 
n(x, y, z) is now expressed as a superposition of modes of the ideal 

waveguide. 

E t = £ a,£„ (14) 

¥ 

and 

Ki-L^Krf ( 15 ) 

The summation symbol in (14) and (15) indicates summation over 
guided modes and integration over radiation modes. Indicating the 
index v by the symbol p in case that it belongs to the continuum of 
radiation modes, we have to replace 

£->x:/"V (is) 

y J 

The sum in front of the integral on the right-hand side of (16) indicates 
a summation over the various types of radiation modes. 

For the derivation of coupled differential equations for the expansion 
coefficients a v and b v we need the orthogonality relations of the modes 
of the ideal waveguide. 

C° /" " e,- (£„ X Xlddxdy = 2^ PS,,. (17) 

J -oo J -» I Mm I 

The asterisks indicate complex conjugation. The symbol 5,„ indicates 
Kronecker's delta for discrete values of v and /x; it is zero if one of the 
indices labels a guided mode while the other labels a radiation mode, 
and it becomes Dirac's delta function if both indices label radiation 
modes. 

The series expansions (14) and (15) are now substituted into the 
equations (6) through (9). Making use of the fact that the mode fields 
satisfy the equations (10) through (13), we obtain 

£ U^ + #A,\ (e, X 3C„ t ) - *'o>e (n 2 - n 2 )a y 8„\= (18) 
and 

? {(t + «*-) (e - x *"> + st "*' 

X [(^-4) (V ' XK " ) ]} = ' (19) 
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We take the scalar product of (18) with fij, and of (19) with 3C^. After 
integration over the infinite cross section, we obtain with the help of 
the orthogonality relation (17) 

-£ + #„a„ = 2I^» a, (20) 



and 



^ + #A = 2 L £„A (21) 



with the coupling coefficients 

R » = w ^r /-I /_1 (n2 ~ n2) 8i ■ Zvtdxdy (22) 



and 

-1 

4*W P |8; 



y — go J —00 



X [(n* ~ 4) V ' X K "] dxdy - (23) 



Equation (23) can be brought into a simpler form with the help of (12) 
and by performing a partial integration 

*" = SF l ~W /I /_1 ^ (n2 - w?) 8 « ' ""** (24) 

Finally, we introduce the amplitudes c^ +) and c£ -) of forward and back- 
ward traveling waves by means of the transformation 

a„ - cl +) e-^ z + cJ _ Vfc« (25) 

and 

6, = cJ +) e-"* J - c< _) e^"». (26) 

Substitution into (20) and (21), addition and subtraction of the result- 
ing equations, and regrouping of terms results in the desired coupled 
wave equations 






= L {XJr' +) cJ +) fl-«*+fr>i + ^-•- ) c<- ) e-*^-^)^}. (28) 



The coupling coefficients are defined as 



+ *^ ■§£•«..}** (29) 
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The factors and superscripts p and q indicate the symbols (+) or (-) 
or the corresponding factors +1 and — 1. The propagation constants 
jS M are positive quantities. The coupled wave equations (27) and (28) 
provide an exact description of the imperfect waveguide in terms of 
normal modes of the perfect guide. The use of normal modes of the 
perfect guide results in the simple general form (29) for the coupling 
coefficients. 

Our coupled mode theory can be applied to any type of waveguide 
problem such as waveguides with refractive index inhomogeneities, 
tapers, or bends. It may be that the expansion in terms of ideal modes 
of the waveguide does not provide the most convenient basis for some 
problems. For the description of tapers, for example, we face the follow- 
ing situation. Consider a waveguide which is perfectly straight and uni- 
form up to a point where its cross section begins to increase. After 
some distance, the taper connects to a uniform waveguide of constant 
cross section. If we use the modes of the smaller guide for our mode 
expansion, we see immediately that the coupling coefficients have non- 
zero values not only on the taper itself but also throughout the entire 
waveguide of larger cross section. Even though our description is 
precise and yields the right answers, it is inconvenient for the problem 
at hand. It would be far better to use so-called local normal modes that 
do not themselves describe wave forms in any real waveguide but corre- 
spond at each point z along the nonuniform guide to the modes of 
a hypothetical uniform guide whose cross section coincides locally with 
that of the waveguide under study. Using local normal modes results 
in coupled wave equations of the form (27) and (28) but with different 
coupling coefficients. In case of the taper, these coupling coefficients 
would be nonzero only on the taper itself but would vanish on the uni- 
form waveguide sections. Local normal modes are obviously better 
suited for the description of tapers. For the description of waves in 
bent waveguides it would be most convenient to use modes that locally 
correspond to a straight waveguide whose axis is tangential to that of 
the actual guide. In addition to these problems of convenience, there 
exist problems of convergence of the series expansions (14) and (15). 
A series expansion in terms of ideal modes may converge more slowly 
than an expansion in terms of local normal modes. However, we shall 
see that we can use the series expansion in terms of ideal modes to treat 
most problems of waveguides with only slight refractive index differ- 
ences between the core and cladding materials. In addition, it is usually 
possible to guess the form of the coupling coefficients of a particular 
expansion from the coupling coefficients for the ideal mode expansion. 
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For the purposes of this paper, we restrict ourselves to the descrip- 
tion of waveguides with piecewise constant refractive index distribu- 
tions and allow only deformations of the cross section of the waveguide 
core. Figure 1 shows a typical waveguide imperfection. The refractive 
index distributions n and n a coincide inside of the core and in the clad- 
ding region. They differ only near the core-cladding boundary. If the 
boundary has moved outward from its ideal position, the index differ- 
ence n 2 — n 2 equals n\ — n\ in the region where the actual core over- 
laps the ideal cladding; it vanishes everywhere else. If the core bound- 
ary has moved inwards, we have n 2 — n\ = — (n\ — n\) in the region 
where the actual cladding overlaps the ideal core and zero values 
everywhere else. The field components that multiply the refractive 
index difference term in (29) are either continuous, if they are tangen- 
tial to the boundary of the ideal waveguide, or they jump by a factor 
(ni/n 2 ) 2 if they are normal to the ideal core boundary. If we restrict 
the discussion to core boundary displacements that are so slight that 
the fields can be considered constant over the region of the displace- 
ment and to weakly guiding fibers with (ni/n 2 — 1) « 1, we obtain 
from (29) 

y- >^ <, />y,»)-«] 

X [p£l- £w + gSl- £„]}d0. (30) 

It is noteworthy that the approximate coupling coefficient (30) is 
identical to the coupling coefficient (13) of Ref. 12 for the local normal 
mode expansion. The only difference in the appearance of these two 
coupling coefficients consists in the fact that the derivative of the 
boundary function instead of the function itself appears in Ref. 12 and 
that the entire expression is divided by /3„ — /3„. It has been explained 
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Fig. 1 — Sketch of a fiber with distorted core-cladding interface. 
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in a different place 13 that it is the Fourier component of the boundary 
function at the spatial frequency 0, — 0„ that determines the coupling 
behavior of the modes. This fact makes it clear that (we write / instead 

of r — a for simplicity) 

1 df (31) 



i(J3, ~ 0J & 



is fully equivalent to the function / itself as far as its effect on mode 
coupling is concerned. We can carry the argument one step further and 
replace / with 

1 d *f (32) 



(0, - 0,) 2 dz> 

If we replace / = r - a in (30) with (32), we obtain a coupling coeffi- 
cient that vanishes in straight, uniform waveguide sections. It is not 
hard to guess that a coupling coefficient of this type belongs to an ex- 
pansion in terms of local normal modes of a hypothetical waveguide 
that is tangential to the curved axis of the actual guide. The modifica- 
tion of (30) that is indicated by (32) is thus particularly suitable for the 
description of mode coupling caused by bends of the waveguide axis. 
A description in terms of ideal modes or even in terms of local normal 
modes of a hypothetical guide with straight axis is unsuitable for a de- 
scription of waveguide bends since it leads to coupling coefficients that 
do not vanish on the straight waveguide section behind the bend. This 
brief discussion shows that it is not hard to modify the coupling coeffi- 
cients of the ideal mode coupling theory to extend it to the case of local 
normal mode expansions of different types. 

III. SIMPLIFIED DESCRIPTION OF GUIDED MODES OF THE FIBER 

A. W. Snyder 8 realized that the modes of round fibers and their 
eigenvalue equations simplify considerably if use is made of the fact 
that (»i/ft 2 — 1) « 1 applies to most fibers of practical interest. D. 
Gloge 9 went one step further and showed that the mode fields become 
simple in appearance if they are expressed in Cartesian instead of 
the more conventional description in cylindrical coordinates. Gloge's 
technique is useful for even more complicated waveguide structures 
such as tubes. 14 

We write down the field expressions for the guided modes of the 
round fiber without derivation. 9 The mode fields can be polarized in 
two mutually orthogonal directions. We have for one polarization in 
the core region for r < a 
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8>> = 






Jy+i(icr) 



sin {v + l)<t> 
— cos (v + l)<f> 



+ /*-i(«cr) 



[ sin ^0 J 



-COS (ir — 1)0/ J Wd ' 

(34) 



X '~ 2F 



Mo . 



/f+i(kt) 



COS 

sin 



[v + 1)0 

[v + 1)0 



3C X = — nA 



< o T , v f COS P0 

— J„((£r) •< . 

[sin i/0 



(36) 



The field in the (infinite) cladding region f or r > a is obtained from the 
field expressions (33) through (36) by replacing the amplitude constant 
A with [J,(Ka)/Hp l) (iya)^A. In addition, we replace k with iy and the 
Bessel function J v {nr) with the Hankel function of the first kind 
Hj, l) (iyr). The parameters k and y are defined as (k 2 = a> 2 e 0J u ) 

k = (n?fc 2 - p 2 )* (37) 

and 

y - (0* - nlk*)i. (38) 

The remaining field components vanish, S x = and K, y = 0. The two 
sets of circular functions that are shown in the field equations are neces- 
sary to obtain a complete set of orthogonal modes. The functions in the 
upper as well as those in the lower position belong together. We have 
used n to indicate n ~ n\ « n 2 . 

The set of guided modes is still not complete unless we also include 
the orthogonal polarization. We have again for r < a(8 y = 3C X = 0) 



6. = 



\Ak 



./, 



+i 



fcos(, + 1)0 

"' [sin (v+ 1)0 



- J-iW 



cos (v — 1)0 

sin (v — 1)0 



S x = AJ v { K r) 



3C Z = 



iAk 
~2h 



cos v<t> 
sin v<f> 



J v+1 (i<r) 



sin (v + 1)0 
— cos (v + 1)0 



+ Jv-iinr) 



3Cu = nA 



[0. 



i e G , , N (cos i/0 
Vo [sin ^0 



sin ( v - 1)011 

— COS (v — 1)0J J 



(39) 
(40) 

(41) 
(42) 
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The field in the cladding is again obtained by the replacements men- 
tioned above. 

The mode amplitudes must be related to the power P (that is the 
same for all the modes). We have for the fields (33) through (36) and 
for the orthogonally polarized field (39) through (42) 



A = 

with 



*L- 7 «P 



e,ira 2 n(nl — nl)k 2 \J v -i{ko)J v+ i{Ka) 
2 for v = 



1 1 for v ^ 



(43) 
(44) 



The eigenvalue equation of these simplified fiber modes is 

J v {ko) H?'(iya) 

With the help of the functional relations of the cylinder functions, it 
is easy to show that (45) is also valid if v — 1 is replaced by v + 1. 

The simplified guided modes listed here are not the same as the usual 
HE and EH modes of the round fiber. It can be shown 915 that the simpli- 
fied modes listed in this paper result from the usual fiber modes as super- 
positions of an HE and an EH mode. The HE,+i.„ and EH,_i, M modes 
have very nearly the same propagation constant, they are almost de- 
generate. Since this degeneracy is not perfect, our simplified modes are 
not modes in the true sense of the word. They decompose into the 
HE v +i,^ and EH„_i lM modes of the round fiber as they travel along the 
waveguide thus changing their shape. A true mode is defined by the 
fact that only its phase changes (in the lossless case) as it travels down 
the guide. However, the approximate modes do form a complete 
orthogonal set of modes and can thus be used to express any field that 
can exist in the fiber. Even after one of the approximate modes has 
decomposed into HE and EH modes, it can again be expressed in terms 
of approximate modes at this point. The fact that the approximate 
modes are not true modes in the usual sense does not limit their useful- 
ness for studying mode conversion and radiation problems. 

The important fundamental HE n mode of the fiber corresponds to 
the lowest-order approximate mode with v = 0. This is a true mode 
that does not decompose as it travels along the waveguide. 

IV. SIMPLIFIED RADIATION MODES OF THE FIBER 

The radiation modes of the fiber can again be simplified by using 
(ni/n 2 — 1) <5C l. 16 There is a slight complication, however. The simpli- 
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fied description of the guided modes was made possible by the fact that 
they are very nearly transverse modes, their transverse components 
being much larger than their longitudinal components. The radiation 
modes are nearly transverse only if their propagation constants are 
nearly = n 2 k. Since the continuous spectrum of radiation modes ex- 
tends from fi = -n 2 k to 0. = n 2 k, only the modes in the immediate 
vicinity of the two end points of this interval are also nearly transverse 
modes. Throughout most of the spectral region of values, the ap- 
proximation corresponding to that for the guided modes does not 
work. However, we can still use the fact that the refractive indices of 
core and cladding are nearly identical and use the radiation modes of 
free space in the region where our mode approximation technique fails. 
A simplified treatment of all the radiation modes is thus also possible. 
The two approximations complement each other. In the region near 
= ±n 2 k, where we use the approximate radiation modes of the guide, 
the free-space radiation modes do not work very well because reflec- 
tions at the core-cladding interface at grazing angles are important. 
Inside of the range, where the waveguide mode approximation 
method fails, we can use the free-space radiation modes with confidence 
since the interface does not cause much reflection for waves passing 
through it at reasonably steep angles. 

We begin by listing the approximate radiation modes of the fiber for 
near ±n 2 k. The field equations are very similar to those of the guided 
modes. In the fiber core at f < a we have 

*--«■•» ft 5} w 






J ' +lM U („ + 1)0 



-'•^{zfciS}] M 

■ e„ , . , [cos v<f>] fiM 

\—J v {ar) \ *\. (49) 

Ho (sin i>(j>\ v ' 

The remaining field components vanish. The propagation constant is 
a continuous variable for radiation modes unrestricted by an eigenvalue 
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equation. The parameter a is defined as 

a = (n\k 2 - /3 2 )* 



(50) 



Instead of specifying the modifications that are required to transform 
the expression of the field inside of the core into the expression for the 
cladding field, we state the field in the region r > a in detail. 



S,= 



iCp 

2/3 



[(ff#i(pr) + DH%i(pr)) I 



S v = C(H?\ P r)+DH?\ P r))l™l* 



3C« — 



iCp 

2k 



e 



sin (v + 1)0 
cos (v + l)<p 

sin (v — 1)0 
— cos {y — 1)0 



cos (v + 1)0 
sin {v + 1)0 



}] 



(51) 
(52) 

(53) 

(54) 



HP and H™ are the Hankel functions of the first and second kind. The 
parameter p is defined as 

p = ( n 2/c 2 - /3 2 )*. (55) 

The amplitude coefficients are 



and 



C = l -^\_aJ ¥+ ,{aa)H?Kpa) - pJ v {<ja)H?Upa)^B 

_ ^jgjpj - pJ v (aa)Hi l U(pa) 
D ~ aJ v+1 {oa)H?\pa) - pM<xa)m 2 Upa) 



(56) 



(57) 



For the field with the orthogonal polarization, we simply state the field 
expressions inside of the core. It should be apparent from inspection of 
(46) through (57) how the field expression in the cladding is obtained 
from that of the core. The relations between the amplitude coefficient 
are the same in either case. We have for r < a 



p %B*\ rcos(v + l)0\ 



[sin (v - 1)0 J J 
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o r> t i \ f COS v6 1 

«.-*/.(") {^4} (59) 

«-t>6{'-«w[-iKS] 

It remains to relate the amplitude coefficient to the power P. Because 
of the continuous mode spectrum, it is not possible to normalize the 
radiation modes with a finite amount of power. The parameter P is 
defined by the relation 

\ f " f " (fip X 3Cj) -ejxdy = £r P8( P - p'). (62) 

The amplitude coefficient 5 is thus [e„ is defined by (36)] 

(^WpP)i 

^aT*\<rJ,,- 1 (<Ta)Hl 1) (pa) - P J,((ra)^ 1 2 1 ( P a) | ' 

It is important to remember that the radiation modes listed so far are 
valid only in the immediate vicinity of = ±n 2 fc. Inside the & range, 
we use the radiation modes of homogeneous space with refractive 
index n 2 . These modes can be expressed in a number of different ways. 
The simplest expressions would result from a plane-wave representa- 
tion. However, for our present purposes, it seems advisable to use field 
expressions that resemble most closely the radiation modes (46) 
through (54) and (58) through (61) in order to achieve continuity of the 
field expressions throughout the entire range. The modes of the 
homogeneous medium (in vacuum we would say free-space modes) are 
simpler than the radiation modes of the fiber, since one expression ap- 
plies throughout all of space. There is no need to treat the fields inside 
and outside of the core separately. The field expressions that satisfy 
our requirements are [p is defined by (55), we use n = n 2 ] 

c ._ iCp [ f sin (v + 1)0 1 

-W,_ l( pr){ !*<*-}>* 11 (64) 

[_ —COS (v — 1)0 J _ 

S* = (65) 
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COS V(p 



S v = CJ v (pr) 



sin vtf> 



_ _ iCp [Jo 
*" ~ 2fc Vmo L 



Jv+i(pr) 



cos (v + 1)0 

sin (v + 1)0 



- J v -i(pr) 



cos (y — 1)0 



(66) 



(67) 



sin (j> — 1)0 

C ^L/(3 , n 2 i;\., . fcosj>0\ 

+'-"{:«: 8})] (68) 

K » - 4>fcU _ T)L /rtlW I "cos (, + 2)4 
The orthogonally polarized modes are 



= iCp\ 

~ 2/3 L 



S a = CJ f frr) 

Sy = 

v _iCp [Jo 



J,+i(pr) 



cos (c + 1)0 

sin {v + 1)0 



cos v<p 

sin »»0 



Jw+i(pr) 



sin (v + 1)0 
— cos (v + 1)0 



(71) 
(72) 

| (73) 



+ J t -i(pr) 



sin (v — 1)0 
— cos (y — 1)0 



C [Z/0 n»*\[, ,,[ sin {v + 2)01 
*• " -iVSU " T)r W l "cos (* + 2)0 j 



3Cu — ~~ 7 






+ 



/j8 n%\( , , . fcos(i> + 2)0 



+ 2)0 



+'-*"{«:£})]• (75) 
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The amplitude coefficient C is related to P 



C = - 



4 V- kpfiP 



(76) 



_ e v ir{^ + n 2 fc 2 ) 

The modes of the homogeneous medium, like all the other modes, are 
mutually orthogonal among each other. With the help of the relations 

Ux)m i U(x) - J, +l (x)H?>(x) = 4- (77a) 

Zirx 

and 

J,(x)H% x {x) - J y+1 (x)H™(x) = - ~ (77b) 

ITX 

it is easy to show that the radiation modes of the homogeneous medium 
and the radiation modes of the round fiber reduce to the same expres- 
sions in the limit | p\ = n 2 k, ni = n 2 = n. 

V. COUPLING COEFFICIENTS FOR CORE-CLADDING IMPERFECTIONS 

We consider a fiber whose core-cladding interface is described by the 
function 

r(x, y,z) = a + f(z) cos (w0 + ^). (78) 

If we choose a different function f(z) and different phase \p for each 
integer m and sum over the second term on the right side, we generate 
a Fourier series which allows us to describe core-cladding imperfections 
of the most general kind. 

We assume that a given mode labeled v is traveling in the waveguide 
and ask for the coupling from this mode to all other guided and radia- 
tion modes. The function f(z) can be separated out from the coupling 
coefficient by defining 

K,, = R„,f(z). (79) 

We have mentioned earlier that the longitudinal field components 
of the guided modes are much smaller than their transverse compo- 
nents. The same statement is true for the radiation modes only if their 
propagation constant /3 is very nearly equal to ±n 2 fc. For |/3| values 
much smaller than n 2 k, the longitudinal field components of the radia- 
tion modes can be as large as or larger than the transverse components. 
The coupling coefficients contain scalar products of the two fields that 
are coupled together. Coupling coefficients that involve at least one 
guided mode are thus determined primarily by the transverse compo- 
nent of both fields, since the product of the longitudinal components is 
small over most of the range of (3 values. We thus neglect the contribu- 
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tion of the longitudinal components and gain the advantage of much 
simpler expressions for the coupling coefficients. The only region of the 
range where the longitudinal components could make a significant 
contribution to the coupling process is near /3 = 0. Outside of a small 
region near /S = 0, our approximation is reasonably accurate. To this 
approximation, modes with orthogonal (transverse) polarization are 
not coupled by core-cladding imperfections of the form (78). 

For coupling between two guided modes v and n we obtain from (30) 
and (79) with the help of the field expressions 

Ve^ 2iank\_ \ ,/,,_i(/c,,a)./,+i(/ca) J M _i(/c M o)^ M+ i(K M a) | ]* 

The factor e„, m can assume the values 4, 2, 1, or 0. It is zero unless 
n = v ± m. Table I shows the values of e» vm for all possible cases. The 
factor e v is defined by (44), k, and ?„ are determined by (37) and (38). 
Coupling between the guided mode v and radiation modes n must be 
described by two different coupling coefficients depending on the value 
of the propagation constant fi of the radiation mode. For values of |/3| 
close to n 2 k, we use the radiation modes of the fiber and obtain 



Zip* f 

Ve,e H 



pi — - I J 7„Vp./,0c,a),/ M (o-a) 

X iTra\J y - 1 (K v a)J v+ i(K,,a)\*\<Tj ll -i(<Ta)H lt W(pa) ~ pJ A^W^pa) | ' 

(81) 

For p values inside the range — n 2 k < fi < n 2 k excluding the end points, 
we use the radiation modes of homogeneous space and find 

p \n { — — 1 ) fcp/3 yJ,(Kji)JApa) 
£( P . fl) = e^m L \ w 2 / J (82) 



kim 



fe t[2(/3 2 + n 2 A; 2 ) \Jr-tMJtiM | ]* 

Use of (77) allows us again to see that (81) and (82) become identical 
in the limit n\ = n 2 = n, |/3 1 = n 2 k. 

The coupling coefficients for the approximate guided and radiation 
modes of the round optical fiber allow us to solve a large number of 
problems involving fibers with core-cladding interface irregularities. 



VI. FAR-FROM-CUTOFF APPROXIMATIONS 

For purposes of multimode operation it is often desirable to have 
simple expressions for the coupling coefficients which are valid far from 
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Table I — Tabulation of the Factor e M , m for All Possible 

Combinations of Angular Field Dependence of the Modes 

and the Core-Cladding Interface Distortion 



e U vm = unless specified otherwise 



incident mode 




spurious mode 


distortion 


COS v<$> 




4 
2 


cos n<f> 

v = fi = m = 
f v = 0, ft = m 


cos m<i> 




6pvm — " 


I n = 0, v = m 








1 


(0 < fi = v ± m 








[0 < fi = m — v 




incident mode 




spurious mode 


distortion 


sin»«£ 


[ 





sin fuf> 

v or fi = 


cos m<t> 




e^m = < 


+ 1 


< fi = v ± m 
< ft = m — v 





incident mode 


spurious mode 


distortion 


cos V<f> 


sinM"£ 

r o n = o 

2 v = 0, fi = m 


sinm<£ 


Gpiim = ] 


1 ft = v 4- m 
— 1 < /I = »* — w 
1 < /t = m — v 





incident mode 
cos*^ 



spurious mode 

.sin ti<t> 
e^m = 



distortion 
cos m<f> 



incident mode 
cos v<f> 



spurious mode 
cos fi<t> 
e u ,m = 



distortion 
sin m<j> 



incident mode 

sin ci/> 



spurious mode 

sin fuf> 
ea»m = 



distortion 
sinmtf. 



cutoff. We obtain such approximations by using the approximation for 
large argument of the Hankel function (7a X> 1) 



\nrya 



rya 
With (83), we obtain from (45) for ya » 1 

yJ v (Ka) = —kJ v -i{ko). 



(83) 



(84) 



We remarked earlier that (45) is also valid if v — 1 is replaced with 
v + 1. We thus also have 



yj f (na) = kJ v +i(kO). 



(85) 
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Multiplying (84) with (85) and taking the square root results in 

yJ v (ica) = k[— J y -i(Ka)J v+ i(Ka)2 i 

= «[|J^i(«o)/h-i(«*)|;P. (86) 

Far from cutoff, (86) allows us to write the coupling coefficient (80) 
between guided modes in the simple form 

£«m> = E5^_5s*_. (87) 

Ve„e M 2iank 

Very far from cutoff, ya — * °o , we obtain the approximate eigenvalue 
equation 

J v {k v o) = (88) 

from (84) or (85). The roots of the Bessel function J v (k v o) thus deter- 
mine the values of n v a that appear in (87). For higher-order modes K„a 
becomes large so that we can approximate 



J > M ~ ^a C0S [ K > a ~ (> + l)l] 



(89) 



Equation (88) requires that the argument of the cosine function equal 
(2N + l)ir/2. This leads to a direct determination of 



K v a 



(» + w + 1) I (90) 



with N = 0, 1, 2 • • • . The equations (87) and (90) provide us with an 
approximate determination of the coupling coefficient between two 
guided modes without the need for solving a transcendental eigen- 
value equation. In a strict sense, we would have to label k and the 
coupling coefficient with N as well as v. We refrain from burdening the 
symbols with too many indices. 

The coupling coefficients between guided and radiation modes can 
similarly be simplified far from cutoff of the guided modes. The far- 
from-cutoff approximation of (81) is 



V ( — — 1 ) ^Kj^iaa) 



KJV.0 = 1^ ^2 / . (Q1) 

Ve»6 M tTa|tt/'p-i(<ra)Hji 1) (pa) — pJ M (o-a)^ 1 i.i(pa) | 

This equation is valid only for | /3 1 values very close to riik. The coeffi- 
cient that describes coupling between guided modes and the radiation 
modes of homogeneous space, eq. (82), leads to the far-from-cutoff 
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approximation 

V \n ( ^ - 1 ) k P Wj^pa) 

R$* = 1^ L XU2 1 J (92) 

Ve^ i[2(0 2 + w 2 fc 2 )]' 

This expression is valid inside the range — n 2 k < /3 < n 2 k but not near 
|/3 1 = n 2 k. For most scattering problems, it is sufficient to use the 
coupling coefficients to the radiation modes of the homogeneous 
medium. Only if the scattering is sharply forward or backward directed 
do we have to use the slightly more complicated expression (91) of the 
coupling coefficients to the true radiation modes of the round fiber. 
In the remainder of the paper, we apply our results to special cases. 

VII. COUPLING CAUSED BY WAVEGUIDE BENDS 

We consider the case of a straight fiber that is connected to a fiber 
section that is bent with a constant radius of curvature and finally con- 
tinues in a straight section. If the curved piece of waveguide causes 
considerable mode conversion, the system of coupled equations (27) 
and (28) must be solved. However, for slight mode conversion, we can 
use the approximation that the incident mode does not change very 
much while power builds up in some of the spurious modes. In this 
case, we obtain the following approximate solution from (27) : 

Cll (L) = c,(0) f L K llv (z)e i ^"-^'dz. (93) 

J 

We assume z = at the beginning of the curved section of length L. 
The description (32) is most appropriate in this case. The second de- 
rivative assumes the constant value \/R, with R being the radius of 
curvature of the circular bend. We thus obtain from (79) [with / 
replaced by (32)] and (93) 



c»(L) 



c(0) 



^ ( S si *- w T (94) 



We obtain R^ from (80) or in its "far-from-cutoff" approximation 
from (87). The integer m appearing in (78) must be set m = 1 in this 
case, since we want to describe a continuous offset of the fiber which 
results in a bend. It is apparent that the amount of power transfer be- 
tween the incident and the spurious mode depends critically on the 
separation between the two propagation constants. The sine factor in 
(94) describes the phasing between the two modes. If the incident and 
spurious modes travel with equal phase velocity, power would be 
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transferred from the incident mode to the spurious mode in proper 
phase so that all the power could be exchanged between the modes. If 
both modes have different phase velocities (and our formula holds only 
in this case), the two modes get out of step so that after some distance 
the power that is fed from the incident mode to the spurious mode tends 
to interfere destructively and destroys the power that has already been 
transferred. The power in the spurious mode thus builds up and decays. 
This process does not involve reconversion of power from the spurious 
to the incident mode, since this process is not included in our perturba- 
tion solution. 

It can be shown that a guided mode of the type (34) produces a ring- 
shaped far-field pattern if it is allowed to radiate out of the end of the 
fiber. The maximum of the ring as seen from the end of the fiber ap- 
pears at an angle 



9,-f (95) 

The circular far-field pattern on a screen is broken up into 2v bright 
dots corresponding to the angular intensity maxima of the field dis- 
tribution in the fiber. The angle 0„ is useful to distinguish guided modes 
experimentally. It may thus be of interest to express (94) in terms of 
this mode angle. Using (87), (95), and 

*"""* -3& -"■*-£* (96) 

we can write (94) in the following form: 

For all practical applications, the separation between the angles 0„ and 
0„ is so small that we can replace them with one angle 0. According to 
(90) we get for the difference 

«,-l,-Af-:fc£ (98) 

if p = v ± 1. We restrict the discussion to coupling between modes 
with the same value of N but adjacent v values. If the results of this 
discussion are implemented, (97) assumes the form 



ejL.1 2 6 fc 2 a 4 w 4 . . , _ _ v L , QQ , 

= -rr _ fiDfifl2 sin 2 (ft - ft) v (99) 



e#? 7r 6 i2 6 2 
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This equation shows that the amount of power transfer caused by 
waveguide bends decreases with increasing mode angle. 

Next we consider random bends. It has been shown in Refs. 3 and 17 
that the exchange of power among randomly coupled modes can be 
described by coupled power equations. The power coupling coefficient 
is given by 

fc„, = |£ M ,| 2 F(/3„ -0„). (100) 

F(/3 M — (3p) is the power spectrum of the coupling function f(z) or of its 
equivalents (31) or (32). Using (32) and writing \/R instead of d 2 f/dz 2 , 
we obtain for the power coupling coefficient for random bends 



h " M.-W/W 



with the power spectrum of the curvature function 

c (i) = (\r L L L m e ""~ , - Ud '\y (102) 

The symbol < ) indicates an ensemble average. 

The guided modes suffer radiation losses even in uniformly bent 
waveguide sections, R(z) = const. 18 - 19 These curvature losses cannot 
be obtained by perturbation theory and thus are not included in our 
discussion. However, our perturbation theory includes mode conver- 
sion losses between guided modes and radiation losses caused by 
changes in the waveguide curvature. 

In terms of the mode angle 6, (102) can be written in the form 



_ el vl 4n 2 fc 2 a 2 

flaw — 4 t» 



(b)- (103) 



The dependence on the mode angle is contained only in the power 
spectrum of the curvature function. 

With the help of the coupling coefficients (91) and (92), radiation 
losses caused by random bends can be calculated. However, the explicit 
expression will not be given here. 

VIII. MODE CONVERSION AND LOSSES DUE TO DISTORTED CORE-CLADDING 
INTERFACES 

Instead of writing down the general formula for the power coupling 
coefficient between guided modes based on (80), we restrict ourselves 
to the far-from-cutoff approximation. From (87), (95), and (100) we 
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obtain 



with 



*--fe33«fc-« (104) 

F(fi, - fi,) = (\j=j Q L KzW'-M'dzfy (105) 



f(z) is the function that appears in (78). Equation (104) is remarkably 
similar to eq. (37) of Ref. 3 which was derived for the slab waveguide 
model. For comparison of the two equations, it is necessary to re- 
member that we have assumed that the angles 9 y and 0,, are much 
smaller than unity and that yd » 1 according to our far-from-cutoff 
approximation. The difference in the position of n in the two equations 
is attributable to the different definitions of the mode angles. In our 
present discussion, we consider 0„ as the angle of the far-field radiation 
cone outside of the fiber, while this angle was denned as the angle of 
the plane waves inside of the fiber core in Ref. 3. The correspondence 
between the two coupling coefficients requires us to consider the case 
of pure diameter changes, m = 0, and assume that both modes have 
no circumferential variation, v = and n — 0. In this case, we have 
e 2 ooo/(±e eo) = 1 instead of the factor 1/2 appearing in (37) of Ref. 3. 
The difference is accounted for if we remember that the round fiber 
corresponds to a slab in which both interfaces have irregularities which 
are perfectly correlated. The slab waveguide theory of Ref. 3 assumed, 
on the other hand, that the two interfaces had uncorrelated irregulari- 
ties. This comparison shows that the results of the slab waveguide 
theory and the round fiber are in very good agreement. Our present 
formula (104) holds for core-cladding interface irregularities of a much 
more general kind. Not only pure diameter changes but elliptical 
deformations and deformations of even more general shapes are 
included. 

Next, we turn to the problem of radiation losses. The power loss 
coefficient is defined by [compare (9.3-14) and (9.3-42) of Ref. 1] 

/7l2fc I ft I 

\K„\*F(0-p,)^d0 (106) 

-n 2 k P 

with Ry., being the coefficient for coupling between a guided mode v 
and a radiation mode with angular symmetry n and propagation con- 
stant /3. The power spectrum F is defined by (105) with /3„ replaced by 
/3. Not much can be gained from substituting (81) and (82) or their 
far-from-cutoff approximations (91) and (92) into (106). The integral 
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in (106) is hard to solve and useful approximations covering the whole 
range of correlation lengths have not yet been found. However, using, 
for example, (81) and (82) in (106) simplifies the numerical integration 
compared to the problem discussed in Ref. 7. The radiation losses can 
be calculated from (81), (82), and (106) with much simpler computer 
programs and at considerable savings compared to the theory of Ref. 
7. We use an exponential correlation function 



R(U) = (f(z)f(z + u)) = a* exp (- \u\/B) 
and obtain (Ref. 1, p. 371) 



(107) 



F{B - ft) = 



[03-/^) 2 + j- 2 ] 



B 



(108) 



The resulting radiation losses for pure diameter changes, m = 0, are 
plotted for the HE n mode, v = 0, in Fig. 2 as functions of the ratio of 
correlation length B over core radius a for ni/n 2 = 1.01. The curves 
were obtained by numerical integration of (106) with R^ of (81) in 
the range 0.95n 2 k ^ |/3| ^ n 2 k and with R^ of (82) in the range 
— 0.95n 2 fc < /3 < 0.95n 2 A;. For small index differences between core 
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Fig. 2 — Normalized radiation losses of the HE n mode, v = 0, as fimctions of the 
ratio of correlation length B to core radius a for different values of ka = 2ira/\ for 
pure diameter changes, m = 0. ni/n 2 = 1.01. The dotted line results from using 
only "free space" radiation modes. 
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and cladding, the losses calculated from our simplified theory are in 
perfect agreement with the theory of Ref. 7. For larger index ratios, our 
approximate theory begins to fail. For ni/n 2 = 1.43, the error caused 
by our approximation is in the order of 60 percent. The simplification 
gained from using (81), (82), and (106) is apparent by glancing at the 
complex formulas of Ref. 7. 

The dotted line in Fig. 2 was computed by using the radiation modes 
of homogeneous space alone so that (82) instead of a combination of 
(81) and (82) was used in (106). It is apparent that the radiation modes 
of homogeneous space are not suitable to calculate the radiation losses 
for large values of B/a. It was pointed out in Ref. 5 that large B/a 
ratios lead to forward scattering. The radiation makes small angles with 
the core-cladding interface so that reflection at this interface becomes 
important. It is thus necessary to use the radiation modes of the fiber 
for values near n 2 k. 

Our theory allows us to calculate radiation losses for more general 
core-cladding interface distortions. As a second example, we consider 
elliptical deformation, m = 2, and plot the result of the numerical inte- 
gration of (106) in Fig. 3. The power spectrum (108) of the function 
/(«) [defined by (78)] was used again. We also used a combination of 
radiation modes of the fiber and of free-space radiation modes in the 
same way as indicated before. The radiation losses caused by elliptical 
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Fig. 3 — HE a mode radiation losses caused by elliptical deformations, m = 2, of 
the core-cladding boundary. ni/n» = 1.01. 
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core-cladding interface irregularities are smaller than those of pure 
diameter changes. 

Actual numerical values of radiation losses obtained from curves 
like Fig. 2 and Fig. 3 were discussed in previous publications. l - 6 - 7 

IX. CONCLUSIONS 

We have presented a simplified theory of mode coupling in imperfect 
round optical fibers. The simplification was a result of restricting the 
discussion to fibers with small values of ni/n 2 — 1. The simplified 
theory results in much simpler expressions for the guided and radiation 
modes of the fiber and consequently leads to simple expressions for the 
coupling coefficients. For small core-cladding index differences, the 
simplified theory is in excellent agreement with more general theories. 

The principal contribution of this paper is a tabulation of coupling 
coefficients for coupling between guided and radiation modes that are 
necessary for solving mode coupling problems caused by general core- 
cladding interface imperfections. A general coupling theory based on 
the modes of the ideal fiber is also presented. 20 
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